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Abstract 

The Macintosh application StdMap allows easy exploration of many of the phenomena 
of area-preserving mappings. This tutorial explains some of these phenomena and presents 
a number of simple experiments centered on the use of this program. 

Dedicated to the memory of John Greene. 
His ideas are as timeless as his generosity is legendary. 

1 Introduction 

Area-preserving mappings give rise to incredibly rich dynamics and mathematics. They also 
arise in many applications including mechanics, chemistry, celestial dynamics, plasma physics, 
condensed matter physics, and many other fields. 

The "standard map" z — ► f{z) = z' , where z = (x, y) £TxR, is defined by 

x' = x + y - J^sm(2Trx) modi, 
V = V - 2i sm(27rx) . 

It was proposed by Boris Chirikov as a paradigm for resonance phenomena in conservative 
systems; it was also derived independently by Bryan Taylor as a model for the motion of charged 
particles in a strong magnetic field. The standard map is one of a class of maps called "area- 
preserving twist maps;" it also has a number of symmetries, the most important of which is a 
time-reversal symmetry, see 

The map % has a single parameter, k, that represents the strength of the nonlinearity. A 
number of physical systems can be modeled by this map. One is the cyclotron particle accelerator 
as described in my review article [1] ; another is the "kicked rotor" , described in [2] , and third is 
the Frenkel-Kontorova model of condensed matter physics studied extensively by [3] . 

The application "StdMap" (available from <http:/ /amath. colorado.edu/~jdm/stdmap.html>) 
allows you to explore the dynamics of a dozen or so reversible, area-preserving mappings — 
including the standard map — on a Macintosh™ computer. It has an advantage over more gen- 
eral computer algebra tools of being compiled and dedicated to this one task (you can iterate 
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hundreds of thousands of steps per second on any recent computer) and of being interactive. In 
this tutorial, I will describe some of the experiments that can be performed with StdMap and 
demonstrate some of the basic phenomena of conservative dynamical systems. 



2 Starting StdMap 



StdMap is a Macintosh application and runs under Mac OS 10.2 or later. The most recent 
version is always available from <http:/ /amath. colorado.edu/~jdm/stdmap.html>; this tutorial 
will refer to version 4.5. The file on the web site is compressed; when you download it, it may 
automatically decompress, or you can double click on the ".zip" file to create the executable 
application^] 
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Figure 1 : Starting up the application StdMap results in continuous iteration of |T]) for randomly selected initial 
conditions. The text window provides information about the speed of your computer, which is used to select the 
global variable N, the number of iterates per function call. 



Double-click on the application to launch it. After a slight delay, two windows will open, as 
shown in Fig. [T] the left window contains a plot window and the right is a text window. By 
default, the program will start iterating ([TJ for k = 0.971635406 with random initial conditions^] 

What you see on your screen, and in Fig. [T] are many orbits, each given a color from a list 
of 19 colors. Each orbit consists of a randomly selected initial point zq = (xQ,yo) and its next 
./V iterates, 

(xt,Vt) = f(x t -i,y t -i) , l<t<N, 

where N is selected by the program depending upon the speed of your computer^ Some of the 
orbits behave "regularly" filling out topological circles or families of circles and some behave 

1 For the purposes of this introduction, if you have ever used StdMap before, you should delete the preference 
file ~/Library/Preferences/edu. Colorado. stdmap.plist, so that you can start the application in its default state. 

2 For the significance of this parameter value, see 

3 On my 3 GHz machine, N is 91074 to give an elapsed time for N iterations of about 1/60-second. You can 
view or change this using the Change menu and selecting Speedometer. . . . 
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"irregularly" or "chaotically" filling out what seems to be a region of nonzero area. We will 
discuss this more in Sj4] 

The natural phase space of the standard map is the cylinder, Txl, with the angle variable, 
x, having period-one, and the momentum, y, taking any real value. To implement this, the 
"mod" operation for x in ([I]) is defined in the program as 

x mod 1 = x — [x + 1 /2j 

so that x mod 1 G [— |, i). However, the standard map is actually periodic in the momentum 
direction as well, and the phase space can also be thought of as the two-dimensional torus if 
we add a mod 1 operator to the y-equation; indeed, this is the default setting for StdMap. The 
reason why this works is that the orbit starting at (xo,yo) and the orbit starting at (^Oj^o) = 
(xo + m, yo + n) for any (m, n) € Z 2 are related by 

6 = x t 

m = yt + n 

thus the vertical separation by an integer stays constant, and r]t mod 1 = yt. You can turn off 
the periodic modulus for y by selecting the Change menu (see Fig. [2| and from the Clipping 
submenu, switching to y unbounded. 

3 Iterating Step- by- Step 

Iteration of a specific initial condition in the standard map is easy: just click on the point in 
the plot window to start there. If you would like to start with a clean window, select Clear 
from the Edit menu or type §£-Y. Now click on a point near the origin and notice that it almost 
instantaneously fills in a small ellipse about the origin. This is because (0, 0) is a stable fixed 
point for this value of k, see §|4j and Moser's famous twist theorem, part of KAM theory, implies 
that the iterates near an elliptic fixed point will generically lie on invariant circles [H SI [5] . 

You can also show the number of iterates on an orbit by selecting Change — > Show # of 
Iterates, and the position of the cursor in the window with Change — > Show Position (§g-U), 
see Fig. [2} These settings will be remembered the next time you start StdMap. 

Now click on an initial condition near the unstable fixed point at (0.5, 0.0). This trajectory is 
chaotic and fills a large region of phase space; what is most surprising, though, is that this region 
is highly nonuniform, containing many holes. Numerical evidence indicates that the invariant set 
that is densely covered by a chaotic orbit is a fat fractal (that is a fractal with positive measure), 
though this has never been proven j6[ U]. How long do you have to iterate until the trace of 
the orbit settles down and no new pixels are filled? If you have a monitor with many pixels, 
this time can be very long indeed [SJ. If we change the parameter value by selecting Change 
—* Map Parameters. . ., and typing the value 2.0 for k in the Parameter Dialog window, the 
large chaotic region has the form shown in Fig. [3} 

Often iteration in StdMap is too quick to really see what is going on; there are two ways 
to slow it down. One is to change the iteration mode in the Find menu. In particular Find 
— ► Single Step (spacebar) stops the iteration and moves the point forward (now drawn as a 
small square instead of a pixel) only when you hit the spacebar, recall Fig. |2j Let us also change 
the parameter of the map to something smaller so that there is less chaos, choose Change — » 
Map Parameters . . . , and type the value 0.3 for k in the Parameter Dialog window. Note that 
the current value of k is displayed at the bottom of the plot window. Now when you click on an 
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Figure 2: The Change and Find menus for StdMap. In this tutorial, we denote a menu selection using the — > 
symbol, thus Change — > Show Axes indicates that we are selecting the sixth item in the change menu. This item 
is checked since the axes are currently shown. Selecting it will toggle the display of the axes in the plot window. 
Your selections in this menu will be remembered the next time you start StdMap. 



initial point, and repeatedly hit the spacebar to iterate step- by-step, you will see a portrait like 
that shown in Fig. 4| 4 



Iterating one step at a time reinforces the fact that maps are dynamical systems with discrete 
time. Another key feature of the map ([!]) is that the horizontal distance between successive 
iterations grows with the momentum value. Mathematically this is an example of the twist 
condition, 

< 2 > 

For the standard map, r = 1, and so it twists to the right. Perhaps a better way to visualize 
twist is to iterate a curve of initial conditions instead of a single point. You can do this in 
StdMap by selecting Find —* Curve ... or typing §§-J, recall Fig. |2j This will open the curve 
dialog, as shown in Fig. [5j There are five types of curves that you can iterate, and you can 
select one by clicking in one of the boxes. For this demonstration, click on the middle box, which 
selects the line type. Now click and drag a line in the plot window. A small dialog window will 
open near the bottom of the plot when you release the mouse. This shows the number of iterates 
that you have performed, either by clicking the Iterate button or hitting the spacebar or the 
I key. If you start with a vertical linej^] as we did in Fig. [HJ the first iterate will be the line 
y = x — the vertical line twists to the right and becomes a graph over x. This is the geometrical 
meaning of ^ . The second iterate of this line becomes the curve 

(x, y) = (2s — — sin(27rs), s — — sin(27rs)) 

27T 27T 

which is still a graph over x for \k\ < 2.0. When k = 0.3, the curve is no longer a graph on 



4 The second way to slow iteration is to decrease the global variable iV that represents the number of iterates 
done each function call when you are in continuous iteration mode, %%-G. You can do this with Change — > 
Speedometer. 

5 It is a standard Macintosh convention that drawing is constrained if you hold down the shift key. In this 
case, the line will be constrained to be exactly vertical, horizontal or at forty-five degrees. Constrained drawing 
is also an easy way to draw squares and circles. 
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Figure 3: Chaotic orbits appear to densely cover a fat fractal like this set generated by iteration of a single 
initial condition for k — 1.0 (left pane) and k — 2.0 (right pane) 



the sixth iterate — it is triple-valued for some intervals of x, see Fig. |5] Nevertheless, this curve 
still has some aspects of twist: it tilts to the right. The map f n is a composition of twist maps, 
and is called a tilt map jS]. Tilt maps have many of the properties of twist maps, and it has 
even been recently shown that they have variational principles, so that KAM and Aubry-Mather 
theory can be applied [TO] . 

4 The Onset of Chaos 

The standard map is integrable when k = 0. Indeed for this value of k, the momentum is an 
invariant, and all the orbits lie on horizontal curves. More generally, an area-preserving map is 
integrable [JJ if there exists a non-constant function I(x, y) such that 

Iof = I. 

There are several other integrable maps that can be studied in StdMap; for example, the McMil- 
lan map, 

{x',y') = (x,-y-ax-b- ^\ , (3) 

is integrable whenever a = 1 or a = 0, for any b |12j . The second is an example of a generalized 
standard map, 

{x',y') = {x + y + F{x),y + F{x)) , (4) 



where the forcing function is 



^. . 2 ( asin(27T2;) , 

Fix) = arctan . . . (5) 

v ; vr V 1 + «cos(2vrx ' 
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Figure 4: Iterating the standard map one step at a time for k — 0.3. 



This map and generalizations of it were proved to be integrable by Suris P3J. Note that (jjjj) 
limits on the standard map as a — > 0, indeed when \a\ < 1, F has a simple Fourier series: 

, , 2 ^ {-a) n , 
F(x) = - 2^ - — sin(2vrnx) 

71 n=l 71 

whose first term is just the standard map with k = 4a. The invariant functions for these maps 
are more complicated than that for the standard map; however, in each case there exists an 
analytic invariant I(x,y). You can select either of these maps and study their dynamics using 
the Mapping menuj^] 

Now return to the standard map at k = 0. Since the phase space is periodic in x, the orbits 
lie on topological circles, y = c; we call them rotational circles since they encircle the cylinder, 
as opposed the librational circles that lie near the origin for k > 0, recall Fig. [TJ What happens 
to the rotational circles as k increases? This is the subject of KAM theory — the theory proposed 
by Kolmogorov and developed by Arnold and Moser. A nice way to visualize this with StdMap is 
to start a trajectory and increment the parameter slowly. 

First set the parameter & to by typing §g-k and entering in the dialog. Now click in the 
plot window to start an orbit (if you are still in single-step mode, type 3€~G to start continuous 
iteration). The up-arrow and down-arrow keys can be used to increment k. The step size, Ak, is 
set using Change — * Parameter Increment . . . ; by default it is set to 0.05. Repeatedly hitting 
the up-arrow key will restart the orbit at the new parameter value using the current point as 
the initial condition. You should see a sequence of pictures like those overlaid in Fig. |6j 

In the left pane of Fig. |6j the initial orbit is selected with a small value of y; when k = 0.1, 
it is deformed to a rotational invariant circle, but when k = 0.2, it is trapped into the main 
resonance about the fixed point at (0,0). This librational invariant circle evolves slowly until 
k = 2.25, when the orbit is trapped in a chaotic layer about a period-4 saddle. In the right 

6 In StdMap, the Suris map has a second term Tvbsm(2irx) that destroys integrability when 6 7^ 0. Type 
to access it; here X is the option key — hold this key down and then type This map was studied in [14] . 
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Figure 5: You can iterate a curve by opening the curve dialog from the Find menu and clicking on one of the 
five curve types. Click and drag in the plot window to create the curve, then click on the Iterate button for each 
step, and finally Stop to finish. The five curve types are (1) small boxes, (2) rectangles, (3) lines, (4) polygons 
and (5) piecewise linear curves. In the last case each click forms a vertex; to make the last vertex, double-click. 
The left pane shows 9 iterates of the vertical line x — —0.1 when k — 0.3. 



pane, the orbit is started with a larger y value, and it remains a rotational invariant circle until 
k = 0.85, when it is trapped in a chaotic layer surrounding the main resonance. Note that the 
steps that you will see in your experiment depend in detail on the initial condition, the step 
size, and the current position (xt,yt) when k is incremented. Nevertheless, the phase space 
of the standard map consists primarily of regular invariant circles (some rotational and some 
librational) when k < 0.5 and primarily of chaotic orbits when k > 3. 

It is also possible to increment k with the up-arrow key while you are in random iteration 
mode, Find — > Random Initial Conditions . . . , or using £N§§-I, where the symbol -O represents 
the shift key in the menus, recall Fig. [2] 

Where does the chaos come from? A partial answer is: through bifurcations of periodic 
orbits. An orbit {(xt, yt) t G Z} is periodic if there is a q G N such that (x q , y q ) = (xo, yo); the 
least q for which this is true is the period of the orbit. For dynamical systems on the cylinder it 
is convenient to partially classify periodic orbits by their rotation number, oj = |, were p is the 
number of times the orbit encircles the cylinder. This can be easily computed by omitting the 
mod 1 in (jTJ), so that the resulting map is a lift of the standard map to M 2 . For the lift, 

an orbit is (p, g)-periodic if (x q ,y q ) = (xq +p,yo)- 

Periodic orbits can be computed in StdMap by specifying the pair (p,q) G N 2 . To do this 
select Find —* Periodic Orbit. . .. For example selecting (p,q) = (2,5) results in the output 

(2,5) at (-0.5,0.262766939830826) (6) 

when k = 2.0. This particular period-5 orbit is a saddle^ We can determine this by turning on 

7 The initial condition dialog has a pop-up menu labeled symmetry. The default setting for this is minimize, 
and this gives the orbit above. Choosing minimax from this menu will give the elliptic (2, 5) orbit. The symmetries 
of periodic orbits will be discussed in 
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the calculation of the residue from the Change menu. John Greene introduced the residue as a 
convenient stability index for area-preserving maps in 1968 [15]. Let Df denote the Jacobian 
matrix of the map; for Q, 

nff \ /l-A:cos(2vrx) l\ 

Df(x,y) = . (7) 

y —k cos(27txj 1 J 

The linearization of / about a period-ra orbit {zq, zi, ■ ■ ■ , z n -x} is the matrix 

n-l 

M=\{ Df(zt) = Df(z n ^)Df(z n ^ 2 ) . . . Df( Zl )Df(z ) , 
t=o 

and the eigenvalues X± of M are the multipliers that determine the linear stability of the orbit. 
Since / is area-preserving, det Df = 1, so that detM = A+A_ = 1. Thus the eigenvalues of M 
are completely determined by its trace, trM = A + + A_. Greene defined the residue as 

fl=i(2-trM). (8) 

Since the product of the two multipliers is 1, it is easy to see that orbits come in four classes: 

1. hyperbolic saddle when R < with multipliers < A_ < 1 < A + ; 

2. elliptic when < R < 1 with multipliers X± = e ±2muJ on the unit circle; 

3. reflection-hyperbolic saddle when R > 1 with multipliers A + < — 1 < A_ < 0; and 

4. parabolic when R = or 1 with a pair of multipliers A+ = A_ = ±1. 
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In particular ([I]) has two fixed points, the origin and (0.5,0) = (-0.5,0). From 0, the 
residue of a fixed point is simply R = | cos(27rx), and so origin has residue 

k 

#(o,o) = ^ , 

implying that it is elliptic when < k < 4 and reflection- hyperbolic for k > 4. This can 
easily be seen in StdMap by incrementing k (using the up-arrow and down-arrow keys) and 
studying the orbits near the origin: when k reaches 4 the origin undergoes a period-doubling 
bifurcation and becomes unstable. Indeed one can show that there is an infinite sequence of 
such period doublings that branch from this orbit; they accumulate in a universal way that was 
first discovered by Feigenbaum for dissipative maps and by Greene and collaborators for the 
area-preserving case [16] . 

By contrast, the residue of the second fixed point is 

k 

#(0.5,0) - -4 ; 

this always negative for k > 0, so that this point is a saddle with positive multipliers. 

Turning on the residue calculation for the (2,5) orbit at k = 2.0 and finding the orbit 
again, gives R = —11.88, confirming that this orbit is also a regular saddle. Indeed, if you 
start iterating at this orbit, by selecting Find — > Continuously (which automatically fills in the 
initial condition of the last orbit) you will find that the instability leads to numerical error: the 
iteration falls off the periodic orbit and rapidly covers the chaotic fat fractal shown in Fig. [3} 
Another way to see this is to find the (2,5) saddle, Find — > Periodic Orbit. . ., and then 
immediately select Find — » Single Step to start single-step iteration at the point ([6]). As you 
iterate with the spacebar, you will see that the orbit visibly deviates from the periodic orbit 
after about 40 steps. For the given residue, the unstable eigenvalue of M is A + = 49.5, so that 

N /5 1 

any error grows by a factor of A, ~ 4 x 10 by N = 40; thus given the inevitable truncation 
error of double-precision floating point computations, it is not surprising that the orbit is lost. 

There are better ways of investigating the properties of hyperbolic saddles, as we describe 
next. 



5 Stable Manifolds 

The stable and unstable sets, W s and W u , of an invariant set A are defined as 

Wl = {(x,y): f(x,y)^ A as i ^ 00} , 
Wt = {(x, y) : f\x, y) -> A as t -> -00} . 

If A is hyperbolic, then these sets are smooth submanifolds that are tangent to the eigenspace 
of the linearization of the map at A |17 [ 118 ) . For example, at the hyperbolic fixed point (0.5, 0), 
the linearization ([7| has eigenvectors 

v± = 

with A± = \ (2 + k ± y/k(k + 4)) . When k > the +- eigenvector corresponds to the unstable 
direction, and since A_ < 1, it has a positive slope. The stable direction has negative slope. The 
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stable manifold theorem implies that W u ' s are smooth curves that start at (0.5, 0) with slopes 

To find these, StdMap uses an algorithm suggested by Dana Hobson |19j . Selecting the 
command Find — > Stable Manifold. . . (^g-H) opens a dialog with three text fields^] The first 
two, p and q, define the rotation number of the periodic orbit that will be the invariant set 
A. The third field is a sign, s = ±1, that determines which branch of the manifold to follow, 
i.e., one of the directions ±v for eigenvector v. The default values (p, q) = (0,1) and s = 1 
will plot the upward-going unstable manifold for the hyperbolic fixed point. After entering the 
three values, the program first finds a point z* on the (p, q) orbit and computes its eigenvectors, 
then it allocates a large block of memory to store the computed points and finally it selects an 
initial point zq = z* + sev+ where e is chosen so that the distance A+e|-u+| is about one pixel. 
Clicking on the Iterate button or hitting spacebar iterates this point q steps. The two points 
{zQ,z q } form an approximation to a fundamental segment of W u for the (p, q) orbit. If these 
points are more than one pixel apart the program inserts a new point between them using linear 
interpolation. 

Hitting spacebar again causes all of the points on the fundamental segment to be iterated q 
more times, and more points to be interpolated as needed. The result is a growing curve forming 
the unstable manifold. The color scheme follows a convention suggested by Bob Easton: red — for 
unstable — represents the blood moving away from the heart, and blue — for stable — represents 
the blood returning. The fundamental segment is stored in the large array and iteration will 
stop if this array fills up. Examples of unstable manifolds for k = 1.0 and 2.0 are shown in Fig. [7} 
Though W u begins life as a nearly straight line along v+, it rather quickly develops intricate 
whorls and tangles, and to the approximation of the computer screen appears to densely cover 
regions of phase space. 

You will notice that stable manifolds (blue curves) are also computed by StdMap — these are 
computed essentially for free using the reversing symmetry, see §|6j What is especially interesting 
is that the stable and unstable manifolds do not coincide, and that when they intersect, they 
usually cross at a nonzero angle; indeed for the (0, 1) orbit this was proven by |20[I21|. Points that 
are on both the stable and unstable manifolds of an invariant set A, were termed homoclinic 
points by Poincare, who discovered them. By the definition Q, homoclinic points are both 
forward and backward asymptotic to A. The complexities of this picture were first envisioned 
by Poincare who famously said: 

When one tries to depict the figure formed by these two curves and their infinity 
of intersections, each of which corresponds to a doubly asymptotic solution, these 
intersections form a kind of net, web or infinitely tight mesh; neither of the two curves 
can ever cross itself, but must fold back on itself in a very complex way in order to 
cross the links of the web infinitely many times. One is struck by the complexity of 
this figure that I am not even attempting to draw |22j . 

You should immediately notice that there is a strong correlation between the chaotic fat 
fractal and the set covered by the invariant manifolds: compare Fig. |7] with Fig. [3| Indeed, 
one reasonable hypothesis is that the closure of W? 5 Q s is identical to the connected chaotic, 
fat fractal containing the saddle fixed point. The transverse intersection of stable and unstable 
manifolds gives rise to the famous horseshoe structure discovered by Stephen Smale in the 1960's 

8 There is also a pop-up menu for selecting the symmetry of the orbit. Leave this set at the default minimizing 
as this corresponds to the hyperbolic saddle. 

9 We allocate an array that is no more than half the physical memory in your computer and no larger than 
needed to contain up to 10 times the number of pixels in the plot window. 
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[23], and the existence of a transverse homoclinic orbit is one of the key features of chaos. As 
an exercise, you can also plot the stable and unstable manifolds for the Suris map for 
the integrable case, b = 0, to see that its stable and unstable manifolds coincide. 10 




Figure 7: Stable and unstable manifolds for the standard map at k = 1.0 (left pane) and k = 2.0 (right pane). 
The left pane shows manifolds for the orbits (1,2), (2,5), (3,8), (1,3), (2,7), (0,1), (-2,7), (-1,3), (-3,8), 
(—2, 5), (— 1, 2), The right pane shows only the manifolds for (0, 1). 



6 Reversing Symmetries 

All of the maps in StdMap are reversible: there is an diffeomorphism S that reverses the map 

foS = Sof- 1 . (10) 

In this section we will explain why reversible maps are easier to study numerically than their 
irreversible cousins. 

Inversion of a map in StdMap is done using (10) rather than writing a separate subroutine 
for / . You can iterate backwards in StdMap by selecting Change — ► Inverse Iteration. 
The inverse will not look very different unless you are doing single steps, using Find — ► Single 
Step, CN3§-G, or iterating curves using Find — ► Curve . . . , g§-J. 

For the generalized standard map Q with an odd force F(—x) = —F(x), one choice for a 
reversing symmetry S is 

Si(x,y) = (-x,y + F(x)) . 
Indeed, it is easy to check that the inverse of Q is 

^{x, y) = Si o f o Si = (x - y,y - F(x - y)) . 



10 Do not set a = 1 for the Suris map, as in this case F(x) = 2x and every point is parabolic and has period-4. 
The Suris map approximates the standard map with a w fc/4, so try a = \ to obtain a picture similar to the left 



pane of Fig. [7] 
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Note that S\o S\ = id, so that S\ is an involution. A consequence is that 

£2 = / Si(x, y) = (y- x, y) 

is also a reversor for /, and is an involution as well. Moreover these reversors are orientation- 
reversing, since their Jacobians have a negative determinant. Thus the generalized standard 
map has been factored into a composition of two orientation- reversing involutions, / = S% o S\. 

One of the significant implications of reversibility is that there are many symmetric periodic 
orbits, namely orbits invariant under S. It is not hard to see that any symmetric orbit must have 
points on the fixed set Fix(S') = {z : S(z) = z}. For the reversors of the generalized standard 
map these are 

Fix(5i) = {(0, 1 /)} > 
Fix(5 2 ) = {(x,y) : x = y} . 

Maps can also have other symmetries that conjugate the map to itself, or commute with the 
map: 

Sof = foS. 

The collection of all reversors and symmetries of a map form a group, called the reversing 
symmetry group [24 . For example, when the force is periodic, F(x + 1) = F(x), then Q 
commutes with the integer rotation^] 

R(x,y) = (x+ l,y) . 

It is easy to see that the composition of any symmetry with a reversor is itself a reversor. Thus 
the maps 

^(x, y) = S 1 o R( x , y) = (-1 - x, y + F(x)) , 
Sa{x, y) = S 2 o R(x, y) = (y - x - 1, y) , 

are also reversors; they give another decomposition, f = S40 S3 The fixed sets of these reversors 
are 

Fix(5 3 ) = {{-\,y)} , 

Fix(5 4 ) = {(x,y) :y = 2x + l} , 

since when the force is both periodic and odd, then F(^) = F(—^) = 0. 

The map / is itself a symmetry, since it commutes with itself. This implies generally that the 
transformations f n o S, where S is a reversor, are also reversors. The set of reversors generated 
by {Si, R, /} form a family. 

You can plot the fixed sets of this family in StdMap with the command Find — ► Symmetry 
Lines . . . , 9g-L. This command will draw the four fixed sets Fix(S'j), and by clicking the Iterate 
button you can iterate generate the fixed sets Fix(/ 2fc o Si): the image of a fixed set is the fixed 
set of another member of the family: 

/(Fix(SO) = {z : f-\z) = S(f-\z))} = Fix(/ 2 o S) . 

11 Here we are really considering the lift of the map to the plane, by undoing the mod 1 operation. It is 
convenient to do this to count the number of rotations of an orbit. Indeed, as was suggested to me by Robert 
MacKay, StdMap uses a structure, ifnumber, for x containing a long and double for the integer and fractional 
parts. 
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Inverse images of the symmetry lines can be found by selecting Change —* Inverse Iteration. 
Several forward and backward iterates of the symmetry lines for ([TJ are shown in Fig. [8j 

The most interesting aspects of the resulting picture are the many intersections of the sym- 
metry sets. It is not hard to see that fixed sets of reversors in a given family intersect at 
symmetric periodic orbits; for example, if z 6 Fix(S') n Fix(/ 2fc o S), then 

z = S(z) = f k S(z) = S(f 2k (S(z))) = r 2k (z) , 

so z has period 2k. Several orbits near such symmetric periodic orbits are shown in Fig. [8} 




Figure 8: Four forward and four backward iterates of the symmetry lines for |TJ at k = 1.3 and several orbits 
near some elliptic, symmetric orbits. 

The generalized standard map also has the inversion 

I{x,y) = (-x,-y) 

as a symmetry when F is odd. This symmetry can be used in StdMap to "mod" the dynamics 
to the upper half-plane, by selecting Change — * Clipping — > y > 0; in this case points with 
y < 0, are reflected to the positive half-plane using /. The inversion symmetry gives rise to two 
additional reversors 

SsO, y) = Si o I(x, y) = (x, -y - F(x)) , 
Sq(x, y) = S 2 o I(x, y) = (x- y, -y) . 

This "extra" symmetry gives rise to some otherwise surprising bifurcations of the librating 
periodic orbits [25] are it has also been profitably used in the study of nontwist maps [26 . 

Finally, the generalized standard map also has a vertical rotation symmetry, V(x, y) = 
(x, y + 1), in the sense that 

foV = RoVof. (12) 

This symmetry can be exploited in StdMap to "mod" orbits to the range — \ < y < \ or in 
conjunction with the inversion I to < y < |. Both of these settings are available on Change 
-> Clipping. 
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Table 1: Symmetry lines containing (p,q), symmetric periodic orbits for a reversible map with 
reversor S\ and discrete rotation symmetry R. Here S2 = f ° S\, S3, = S\ o R and S4 = S2 R- 
The initial point zq on the symmetry line shown maps to the point z^ on the second symmetry 
line after I rotations. 



Since the fixed sets of the orientation-reversing involutions are one-dimensional, it is easier 
to search for symmetric periodic orbits than general orbits. StdMap exploits this by using a 
one-dimensional secant method to find symmetric orbits when you select Find — * Periodic 
Orbit .... To find an ^-symmetric period-^ orbit one could look for fixed point of f q that lies 
on Fix(S'). However, it is a bit better to halve the work, as follows. Suppose z is a (p, q), orbit, 
i.e., after q iterates, the orbit rotates p times so f q {z) = R p {z). If, in addition z £ Fix(S') and q 
is even then 

/«/*(*) = p'\S{z)) = S{r q/ \z)) = S{R-v{P'\z))) 

which implies that Zk £ Fix(S" o R~~ p ), where k = q/2. Thus to search for such an orbit we 
need to iterate only k times starting on Fix(S) and ending on Fix(S' o R~ p ). When q is even p 
should be odd, since (p, q) should have no common factors. Moreover, since Si flips the sign of 
x, Fix(Si o R-p) = J R(P+ 1 )/ 2 Fix(S'i o R), that is, the fixed set Fix(S 3 ) rotated (p + l)/2 times. 
The full list of possible symmetric orbits depending upon whether p and q are odd or even is 
given in Table [TJ For each (p, q) there are two different symmetric orbits, and they each have 
points on two of the symmetry sets. 

When you select Find — ► Periodic Orbit . . . , the dialog box has a symmetry pop-up menu 
where you can select the beginning symmetry line for the (p, q) orbit; this menu has six entries. 
The last four entries allow you to select one of the four fixed sets for the initial point of the orbit. 
A number of symmetric orbits are shown in Fig. |9] As Table [T] shows, one of the symmetric 
orbits always has a point on Fix(Si); it can be shown that this orbit has positive residue when 



k > 0: Fix(Si) is called the dominant symmetry line for the standard map 12 All of the elliptic 
islands for the rotational orbits line-up on this fixed set. The second orbit is always hyperbolic, 
and either has a point on Fix(5i o R) = Fix(5 3 ) when q is odd, or Fix(/ o Si) = Fix(5 2 ) when 
q is even. 



12 When k < the dominant line for |T]) is Fix(S l 3). Other reversible maps, like the Henon map, %%-2, also 
appear to have dominant symmetry lines, though as far as I know this has not been proved. 
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The first two entries in the symmetry menu, Minimizing and Minimax refer to the action- 
minimizing and minimax orbits of Aubry-Mather theory, see [I] for a review. The minimizing 
orbit is always hyperbolic [27] and corresponds to the second orbit above. The minimax orbit 
has positive residue and, when k is positive and small enough, it is elliptic. It has a point on 
the dominant symmetry line. 




-0.5 x 0.5 



Figure 9: Symmetry lines and some symmetric orbits for the standard map with k = 1.0. 



7 The Critical Golden Circle 

The (0, 1) resonance surrounding the elliptic fixed point is enclosed by a connected chaotic region 
generated by the unstable manifold of the (0, 1) saddle orbit. Chirikov noticed that this chaotic 
region appears to be bounded when \k\ < 1 and unbounded for larger values of |fc| |28| . One 
way to see this is to look for climbing orbits, that is orbits that move from y = a to y = b 
for a < b. StdMap provides a convenient way to do this experiment: select Find — > Transit 
Time. . .. The program then asks you to drag the mouse over a rectangle, IZi, in which the 
initial conditions will be selected — for example, choose a rectangle near the fixed point (0.5,0). 
A dialog window will open that shows the boundaries of the rectangle that you selected and also 
has a field labeled N, which will be the maximum length of any attempted transiting orbits; 
by default, N = 10 5 . Clicking on the OK button will accept these values, and you will be asked 
to drag over a rectangle IZf defining the final region. Choose a small rectangle near the point 



on the (1,2) (minimizing) saddle orbit. 13 After you accept the values for the final rectangle, 
StdMap will iterate randomly selected initial conditions in IZi until they happen to fall in IZf, 
or until the number of iterates reaches N. 



Near, e.g., (0.306,0.612) when k = 1.5. 
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The average transit time 
3.7 x 10 4 , and when k = 1.0 



When you click on the Stop button, StdMap will provide you with a binned list of the transit 
times. While StdMap does not make a plot of this histogram (file a bug report!) it is easy to 
copy the data from the text window (Select it and type §§-C). It can be imported into another 
program such as Maple™ or Matlab™ to make the plot, see Fig. 10 
for k = 1.5 is tt ran s = 501 steps, at k = 1.1 is has grown to ttran, 
it has become 3.9 x 10 6 . In each case the transit-time distribution looks roughly exponential, 
but experience shows that the long-time distribution will decay only algebraically with time 
|221 EQl [T]. Since the average transit time grows rapidly as k decreases towards one. 14 this 
technique is not a good one in which to determine the critical parameter value, k cr , at which 
Urans first becomes oo. 




Figure 10: Histogram of transit times for orbits beginning near the (0,1) saddle and ending near the (1,2) 
saddle for k = 1.5 (black) and k — 1.1 (red). 



A much better way was discovered by John Greene |32 | 1331 34J. This is to look for the "last" 
rotational invariant circle that separates the rectangles IZi and IZf. Indeed, Birkhoff showed 
that if there is no orbit that reaches H t then there must indeed be such an invariant circle [1] , 
and Mather provided a variational construction of climbing orbits [35] . 

Greene's insight was to look at the stability properties of periodic orbits that approximate 
an invariant circle. His hypothesis was that when these nearby orbits were "stable", the circle 
should exist, and when they were unstable, it should be destroyed. A persistent invariant circle 
has an irrational rotation number u and, as we learn from KAM theory, uj should satisfy a 
Diophantine condition |T] . Greene used the rotation number to select a family of periodic orbits 
defined as the convergents of the continued fraction expansion of uj. That is, suppose that 



uj = ao + 



1 



a\ + 



[a ,ai,a 2 , ■ . .] 



where at G N are the continued fraction elements. Then 



to. 



Pi 



[ao,ai, . . . ,Oi] 



is the ith convergent to uj. Greene computed the residues Ri of these the (pi,qi) orbits, and 
observed that when k < k cr (uj), Ri — > 0, and the invariant circle exists. When k > k cr (uj), 
Ri —* oo and the circle is destroyed. At the critical parameter, Ri remains bounded and the 
invariant circle exists but appears to be nonsmooth. 



The time grows as tt 



75 {J\> ^cr) 



31 
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At the time, Greene's observation seemed quite preposterous to mathematicians familiar 
with the number theoretic delicacies of KAM theory. Nevertheless, Greene's criterion works 
extremely well and many aspects of his observations were subsequently proven |36| 157] . 

Greene provided strong numerical evidence that the last rotational invariant circle for the 
standard map has to = 7, where 7 = ^ is the golden mean. The symmetries of ([I]), discussed 
in imply that circles with rotation numbers ±u; + m for any m G Z are equivalent. This 
means, for example that the circles 7 -1 = 7 — 1 and 7~ 2 = 2 — 7 are also destroyed at & C r(7)- 

This golden mean is distinguished by its continued fraction expansion 7 = [1, 1, 1, 1, . . .]; it is, 
in the Diophantine sense, the most irrational number. Percival called numbers whose continued 
fractions have an infinite tail of l's, noble numbers, and it appears that the noble invariant 
circles are locally the most robust circles |38[ l39l l4"0] . 

To find a sequence of convergent periodic orbits in StdMap, use the Find —* Farey Path. . . 
command. The dialog window that opens will ask for a pair of numbers (po, qo) and (pi, qi) that 
will form the base of a Farey tree; these numbers must be neighbors in the sense that 

PiQo -PoQi = 1 • 

The default is to choose (0, 1) and (1, 1), which will allow you to construct any number < to < 1. 
We start with the t 
upper neighbor and 



We start with the triplet {u>l, u>c, lor} where u>l = ^ is the lower neighbor, lor = ^ is the 



_ Po + Pi 

LUC = UL®UR= , 

qo + qi 

where © is the Farey sum operation. 

The Farey path of a number u £ (lvl,lor) is the unique sequence of left and right turns that 
lead to to using the recursive algorithm: 

1. If to = cue stop. 

2. If to € (lol,u>c) then the step is L, the next rotation number is lo' c = wl © loq , and the 
triplet becomes {lol, w' c , u>c}- 

3. If to G (u>c lor) then the step is R, lo' c = u>c © lor , and the triplet is {u>Ci w' c , lor}. 

4. Repeat. 

In each case the outer pair of the triplet are Farey neighbors. The path truncates if and only 
if to is rational. It is not hard to see that noble irrationals have a Farey path that eventually 
oscillates . . . LRLRLRLR .... 

Typing a sequence of L's and i?'s (or l's and r's or 0's and l's) in the text box for the Path, 
instructs StdMap to find each of the orbits with rotation numbers on that path. If you turn on 
the computation of the residue, Change — ► Compute Residue, then the Residue is calculated as 
well. The default settings, u>l = j, lor = \ and the path of alternating LK's gives a sequence 
that converges to to = 7 -2 ; the results are shown in Table [2} Here we have computed the 



minimax periodic orbits that sit on the line Fix(S'i), recall (11). 

The values of yt for k cr appear to converge to 0.405086089 as i — > 00. This convergence is 
geometrical and can be understood by a renormalization argument [36]. The residues for k = 
0.971 do appear to converge to 0, though the last value is not believable; careful computations 
show that the minimax orbits always have positive residues. Similarly, the residues for k = 0.972 
grow, and those for k = k cr appear to be bounded, and moreover to converge to 0.25009 |36j : 
the last couple of values in Table [2] again have numerical convergence problems. 
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Table 2: Residues of the convergents of the minimax periodic orbits, with xo = on the Farey 
path for the rotation number 7 -2 for three values of k, including k cr = 0.97163540631. 
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When k > k cr the invariant circle no longer exists, but the beautiful theory of Aubry and 
Mather [3j E] implies that there still is an invariant set with the same rotation number; this 
set is a Cantor set, and Percival called it a cantorus. Cantori can seen in StdMap using the 
Farey path command as above and setting k > k cr . Several examples are shown for the rotation 
number 7 -2 in Fig. 11 The gaps in the Cantor set appear to all be generated by images of 
a single gap. The largest gap is centered on the dominant symmetry line, x = 0, where all 
the elliptic orbits line up. This perhaps visualizes in the most precise way Chirikov's concept 
of resonance overlap [28J: the islands due to the different periodic orbits that approximate the 
irrational rotation number squeeze the invariant circle most strongly where they line up. 

It is also easy to find cantori for other rotation numbers by changing the path in the Farey 
sequence. Since noble numbers are most robust, the cantori with the smallest gaps will be found 
by selecting paths that have a short head sequence, followed by an alternating LR-tail. 



^ ^ ^ ^ ^ ' "' ' 











Figure 11: Golden mean invariant circle and approximations to the cantori for the standard map at k — k c 
0.974, 0.977, 0.980, 0.985, 0.990, 1.0 and 1.05. The cantori have been shifted vertically to avoid overlap. 



When an invariant circle is destroyed, the resulting cantorus has a nonzero flux of crossing 
orbits. This flux can be computed by taking the difference between the action of the minimax and 
minimizing orbits that approximate the invariant Cantor set, see [1] for a discussion. The action 
of orbits can be computed in StdMap by selecting Change — ► Compute Residue & Action. 



8 Exploring other Maps 

StdMap has thirteen reversible, area-preserving maps in the Mapping menu. If your favorite 



system is not included, do let me know and I will consider adding it. 15 Several representative 



phase portraits are shown in Fig. 12 The last two dynamical systems in the Mapping menu are 



actually Poincare sections of Hamiltonian flows. For example, the two-wave Hamiltonian is 

H(x, y, t) = -y 2 k [a cos(27rx) + b cos(2ir(kx — t))] 

2 2tt 

where k G Z. This system has one and a half degrees of freedom, and as such is most properly 
described in a three-dimensional phase space (x,y,t) eTxIx T. The Poincare map (x,y) — * 
(x',y') is the time-one map constructed by integrating the equations of motion for an initial 
condition (x(0), y(0), 0) to t = 1, and plotting the resulting values (x(l),y(l)). In StdMap, this 
integration is done with a symplectic splitting method that is second-order accurate with a time 



15 It would be best, of course, for me to add a parser to StdMap so that arbitrary maps could be iterated. 
However, much of the code depends upon the reversibility assumption, so this would not be easy. 
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step of 0.1. It gives results that are qualitatively good, but should not be trusted for quantitative 
accuracy. I used it to generate the figures in §9.16 of [IS] . 




Figure 12: Four representative phase portraits generated by StdMap. Which of the thirteen systems do they 
represent? 



9 For Help 

StdMap has many other features that are not documented in this tutorial. In many cases the 
menu items should be self-explanatory; if not there are two help facilities. One is the Help menu 
which gives access to several files that briefly describe the features. A second is the Help Tip 
that pops-up if you let the mouse hover over a menu or button in the program. In some cases 
these tips have extended help that can be accessed by pressing the §g key. 

More discussion of the mathematics behind the standard map can be found in [I] 22] [M] H3] 

02]. 
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